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Abstract. A proof that new massive gravity — the massive 3D gravity model 
proposed by Bergshoeff, Hohm and Townsend (BHT) — is the only unitary system 
• at the tree level that can be constructed by augmenting planar gravity through the 

i <~| curvature-squared terms, is presented. Two interesting gravitational properties of the 

BHT model, namely, time dilation and time delay, which have no counterpart in the 
usual Einstein 3D gravity, are analyzed as well. 
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1. Introduction 

For too long physicists believed that gravity models containing fourth- (or higher-) 
derivatives of the metric were doomed to failure by virtue of one detail: they entail 
unphysical ghost states of negative norm. The pure scalar curvature models, i.e., the 
fourth-order gravity systems with Lagrangian C = R + aR 2 and which are tree-level 
unitary, seemed to be the only exception to this rule. Actually, these systems are 
conformally equivalent to Einstein gravity with a scalar field [TJ. Consequently, despite 
having fourth derivatives at the metric level, these models are ultimately second order 
in their scalar-tensor versions. It is, therefore, perfectly understandable that just about 
two years ago the physical community were absolutely amazed to learn that a particular 
higher-derivative extension of 3D general relativity — that is ghost-free at the tree level 
- has been found out by Bergshoeff, Hohm and Townsend (BHT) [2-14]. It was argued 
that this massive 3D gravity model, that is also known as "new massive gravity" , is both 
unitary and power-counting UV finite in its pure quadratic curvature limit |15j, which, 
as it was pointed out by Ahmedov and Aliev [10], violates the standard paradigm of 
its "cousins" in four dimensions [16] . New massive gravity is defined by the Lagrangian 
density 



where k 2 = 32nG, with G being the 3D analog of Newton's constant, and m 2 (> 0) 
is a mass parameter. It is worth noticing that the Lagrangian density given in (1) 
has a reversed Einstein-Hilbert (EH) term. On the other hand, a formal proof of the 
equivalence of the linearized version of the BHT model and the Einstein-Hilbert-Pauli- 
Fierz gravity was given in [2]; incidentally, this proof was reviewed in [T7j. Nevertheless, 
the physical meaning of this equivalence is somewhat unclear; indeed, the linearized 
version of the BHT system is background diffeomorphism invariant, while the Pauli-Fierz 
theory is only invariant under the Killing symmetries of the spacetime (in particular, the 
3D Minkowski space), which clearly shows that a better understanding of the symmetries 
is still lacking [9]. And what about the odd sign change of the EH term previously 
mentioned? At the linearized level, Deser [15] showed that the EH term breaks the Weyl 
invariance of the BHT model without the EH term and, consequently, is responsible for 
giving mass to the graviton. In other words, the higher- derivative terms provide the 
kinetic energy, whereas the EH term provides the mass in this linearized model, thus 
explaining the weird sign change of the EH term. It is remarkable that the EH term 
gives origin to the mass in the linearized version of the BHT system by breaking the 
Weyl invariance and not the expected diffeomophism invariance [9]. 

At this point it would be interesting to ask ourselves about the reason for doing 
research on massive gravitons. The increased interest in recent years in this subject is 
motivated, on the one hand, by the discovery of cosmic acceleration, which might be 
explained in terms of an infrared modification of general relativity that gives the graviton 
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a small mass [IS]; on the other, by the conjecture that some theory involving massive 
gravitons could be the low energy limit of a noncritical string-theory underlying QED 
[19] . As it is often done for so many other gravitational physical issues, it is advisable 
to consider first the possibilities for massive gravitons in the simpler context of a 3D 
spacetime p2]. The BHT model is accordingly the ideal arena for such investigations. 
Our aim in this paper is twofold. 

(i) To show that the BHT gravity is the only tree-level unitary model that can be 
constructed in 3D by judiciously combining the Ricci scalar R with the curvature- 
squared terms R 2 and R 2 . 

(ii) To explore some interesting properties of this remarkable model that have no 
counterpart in the usual Einstein gravity in three dimensional spacetime. 

We describe in the following the steps we shall take in order to accomplish these 
objectives. We start off our analysis by considering in section 2 the most general 
three-dimensional theory obtained by augmenting planar gravity through the curvature- 
squared terms. Now, taking into account that in three dimensions both the curvature 
tensor and the Ricci tensor have the same number of components [20] , we come to the 
conclusion that the Lagrangian density for the theory at hand can be written as 

C = ^(^R + ^R 2 + ^y (2) 

where a is a convenient parameter that can take the values +1 (EH term with the 
standard sign), -1 (EH term with the "wrong sign"), and a and /3 are free coefficients. 
Note that the constants k, a, and j3, have the mass dimension [k] = —\ and 
[a] = [(3} = —1, in fundamental units. We prove afterward that the BHT model is 
the only unitary system at the tree level that can be built from the Lagrangian given in 
Eq. (2). In section 3 it is shown that, unlike what occurs in 3D general relativity, clocks 
are slowed down in a gravitational field described by the BHT model. This gravitational 
time dilation is the basis of the gravitational spectral shift. An expression for a new- 
massive-gravity-induced time delay is obtained in section 4. Finally, we present in 
section 5 some comments and observations. 

We employ natural units, c—H—1, and our Minkowski metric is diag(+l, -1, -1). 
Our Ricci tensor is defined by R^ u = R x = d u T x flX — d\T x flu + ... . A prescription 
for computing the graviton propagator, as well as a list of some identities that greatly 
facilitate this task, are collected in Appendix A. The derivation of an important result for 
checking the tree-level unitarity of a generic 3D gravity model is sketched in Appendix 
B. 

2. Finding a class of tree-level unitary massive 3D gravity models 

To probe the unitarity at the tree level of the models defined by Eq. (2), we make use of 
an uncomplicated and easily handling algorithm that converts the task of checking the 
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unitarity, which is in general a time-consuming work, into a straightforward algebraic 
exercise. The prescription consists basically in saturating the propagator with external 
conserved currents, compatible with the symmetries of the system, and in examining 
afterwards the residues of the saturated propagator (SP) at each simple pole. Let us 
then compute the propagator for the gravity model in Eq. (2). To do that, we recall 
that for small fluctuations around the Minkowski metric 77, the full metric assumes the 
form 

fiv = ?/> + K V ( 3 ) 

Linearizing Eq. (2) via Eq. (3) and adding to the result the gauge-fixing Lagrangian 
density, C g f = ^x(d^ u ) 2 , where 7^ = h^ u — ^n^h, that corresponds to the de Donder 
gauge, we find 

C = jv^, (4) 
where, in momentum space, 



O = 



ok 2 + 



P( 2 ) + ^lp(l) + !f_p{0-w) _ p(O-sw) 

2A 4A 4 A 



_ ^ k2 p(O-ws) , 

4 A 



^--ok 2 + 2an 2 k 4 + -p K 2 k 4 
2A 4 
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(5) 



Here P^ 2 \ pW, p(°-™), p(o-*) > p(o-««) and p{o- ws ) are the ugual three-dimensional 
Barnes- Rivers operators (see Appendix A). 

Therefore, the propagator is given by (see Appendix A) 



O 



-1 



k 2 + 



,p(2) + 



+ 
+ 



y/2 



-ok 2 + 2a/t 2 A; 4 + \K 2 k A fi 
-AAo + 2 + 8Aan 2 k 2 + 3Af3n 2 k 2 



-ok 2 + 2aK 2 k 4 + \n 2 k A P 
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Contracting now the above propagator with conserved currents T tlu (k), {k^ 1 T' lv 
k v T^ = 0), yields 



(6) 
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\t 2 , (7) 



where m 2 , = —j^z, tUq = ^ 8a ^^ K 2 ■ Assuming that there are no tachyons in the model, 
we promptly find the following constraints 



o 



8a + 3/3 



> 0. 



(8) 
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On the other hand, the residues of SP at the poles k 2 = m 2 ,, k 2 = 0, and k 2 
are, respectively, 



m, 



o 



Res(SP) \ k2=m , = - ±{t$, - V) \ k2=m% , (9) 
Res(SP) | fe2=0 =\(Tl v -T 2 ) | fc2=0 , (10) 

itea(SP) | fc2=m2 = i-(T 2 ) | fc2=mg . (11) 

Now, as is well-known, the tree-level unitarity of a generic model is assured if the 
residue at each simple pole of SP is > 0. Keeping in mind that [t 2 u — \T 2 ^j \ k 2 =m 2 > 

and (t 2 v — T 2 ^j \ k 2 =0 = (see Appendix B), we arrive at the conclusion that: 
(i)-Res(SP) \ k 2 =m 2 > if a = — 1 (which implies (3 > and a < 0), and (ii) 
-Res(SP) \ k 2 =0 = 0. Consequently, we need not worry about these poles; the troublesome 
one is k 2 = m\ since i?es(SP) | fc 2 =m 2 < 0. A way out of this difficult it is to consider the 
m — > oo limit of the model under discussion, which leads us to conclude that a = — 
Accordingly, the class of models defined by the Lagrangian density 



2R + P(r 2 V 



(12) 



are ghost-free at the tree level. For the sake of convenience, we replace (3 with 
where m<i is a mass parameter. The resulting Lagrangian density, 



K^rn 
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-R 2 f r2 _ 3 r2 



2 ^2™2 \ "\iv 
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(13) 



is nothing but the BHT model for massive 3D gravity. 

It is worth noting that it is not clear at all whether or not the particular ratio 
between a and f3 we have previously found will survive renormalization at a given loop 
level, even at one-loop; in other words, unitarity beyond tree level has to be checked 
[9]. Most likely the BHT model is nonrenormalizable since it improves only the spin-2 
projections of the propagator but not the spin-0 projection [21J. 

3. Gravitational time dilation 

Einstein 3D gravity is trivial outside the sources; consequently, no gravitational time 
dilation, or slowing down of clocks can take place in its framework. This can easily be 
shown in the particular case of a spherically symmetric distribution of mass M whose 
metric tensor is approximately given by 
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ds 2 = dt 2 - (1 + A) {dr 2 + r 2 dd 2 ) , (15) 

where A = 8GM hi^, with r being an infrared regulator, and r and 9 are the usual 
polar coordinates. 

Introducing now new radial (r') and angular (#') coordinates through the change 
of variables 



(l_A)r 2 = (1-8GM)/ 2 , 9' = (1 - AGM)9, 

we obtain, to linear order in GM, 

ds 2 = dt 2 -dr' 2 -r' 2 d6' 2 . (16) 

The geometry around the spherically symmetric distribution is, therefore, locally 
identical to that of a flat spacetime as it should; however, it is not globally Minkowskian 
since the angle 9' varies in the range < 9' < 2n(l — 4GM). Accordingly, the three- 
dimensional metric (16) describes a conical space with a wedge of angular size equal 
to 8nGM removed and the two faces of the wedge identified. We thus come to the 
conclusion that in the framework of Einstein 3D gravity no gravitational spectral-shift 
occurs due to the presence of the mentioned odd geometrical effect. It is worth noticing 
that in this context, the non existence of a time dilation does not imply that the 
spacetime is necessarily flat; in other words, the time dilation is not a "classical test" 
of 3D general relativity. As we shall see in the following, the aforementioned bizarre 
geometrical effect does not take place in new massive gravity. To do that we have to 
solve beforehand the linearized field equations related to the BHT system. 
The field equations concerning the Lagrangian density 



where Cm is the Lagrangian density for matter, are 



(17) 



2^-%9ttv — ^V^V^i? — 2R^ p \ u R px — -g^DR + DR^ U 



— j^R 2 9nv + -^RR^u 



-T 

4 ^' 



(18) 



where T^ u is the energy-momentum tensor, and G^ v = R^ u — \g^ v R is the Einstein 
tensor. 

The corresponding linearized field equations are given by 



1- 



□ 
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2 » v+ Ak 



k(T 



4 V 2 



V^-T^ ,(19) 
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where rM = k [±Dfc - T ^J, r, = (l + ^)d p7 / + %^, 7„ = V" 1^- Note 
that here indices are raised (lowered) using 77^(77^). 

Mimicking Teyssandier's work on 4D higher-derivative gravity |22j , it can be shown 
that it is always possible to choose a coordinate system such that the gauge conditions, 
T M = 0, on the linearized metric, hold. Assuming that these conditions are satisfied, it 
is straightforward to show that the general solution of (19) is given by 

V = Vv-^ } , (20) 

where h^J is the solution of the linearized Einstein equation in the de Donder gauge, 
i.e., 



nhfJ = ^{Tn» v -T» u l P>yM=0, (21) 
where 7^ = hffi — ^rj^h^, while ip^ v satisfies the equation 

(□ + m*)W = —{T^ - i V T). (22) 

(V) 

It is worth noticing that in this very special gauge the equations for tp^ and h\J are 
totally decoupled. As a result, the general solution to the equation (19) reduces to a 
linear combination of the solutions of the aforementioned equations. 

Solving Eqs. (21) and (22) for a point-like particle of mass M located at r = 0, we 

find 



^00 



kM 



in 



-K (m 2 r) 



h 



22 



~8n 



Ko(m2r) + 2 In 



r 



(23) 
(24) 



where K is the modified Bessel function of order zero. Note that K Q (x) behaves as 
-Inx at the origin and as x~^e~ x asymptotically. Thence, the metric tensor and the 
spacetime interval are given, respectively, by 



/ 1 - 4MGK (m 2 r) 




-[1 + 4GM(K (m 2 r) 
+2 In-!-)] 



\ 










-[1 + 4GM(K (m 2 r) 
+21n^)] 



\ 
/ 



(25) 



ds 2 = [1 - AMGK (m 2 r)]dt 2 - 
+ r 2 d9 2 ). 



1+AGM K (m 2 r) + 2 In 



r 



(dr 2 



(26) 



In the m 2 — > 00 limit, (25) and (26) reproduce (14) and (15), in this order, as 
expected. The geometry around the point-particle is, of course, not locally identical 
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to that of a Minkwoskian spacetime, signaling in this way the possibility of occurrence 
of gravitational spectral shift. Let us then show that the gravitational time dilation 
does occur in the BHT model. 

Suppose that a signal sent from an emitter at a fixed point (r_g, Oe) is received, 
after traveling along a null geodesic, by a receiver at a fixed point (r R , 9 r ) (see Fig. 
1). Now, the difference t R — t E , where is the coordinate time of emission and t R the 
coordinate time of reception, is the same for all signs sent — the wordlines of successive 
signals are nothing but copies of successive signals merely shifted in time. As a result, 
if the t-time difference between a signal and the next is (He at the departure point, the 
corresponding t-time difference at the the arrival point is necessarily the same. However, 
the clock of an observer situated at the point of emission records proper time (r) and 
not coordinate time (t). Accordingly, oIte = a/1 — 4:MGK (m 2 rE)dtE, and similarly 
dr R = a/1 — 4MGK (m 2 r R )dt R . Since (He — dt R , we promptly obtain 



dr R _ y/l-AMGK Q (m 2 r R ) 
dr E y/1 - AMGK (m 2 r E ) 

« 1 - 2MGK (m 2 r R ) + 2MGK (m 2 r E ) 

= 1 + V R - V E , (27) 

where V(r) = %h 00 (r) = —2MGK (m 2 r) is the gravitational potential. This shows that 
if the clock at (r R , 9 R ) is at a lower potential than the clock at (r^, 9e), i-e., V R < V E , 
then dr R is smaller than drE- In other words, the clock that is deeper in the gravitational 
potential runs slower. Eq. (27) is the gravitational time-dilation formula, or redshift 
formula. It is worth noticing that dr R — > drE as m 2 — > oo, implying that no gravitational 
time dilation takes place in the framework of 3D general relativity, which totally agrees 
with the result we have previously found. 

On the other hand, if the emitter is a pulsating atom which in the proper time 
interval Ate emits n pulses, an observer situated at the emitter will assign to the atom 
a frequency ve = which, of course, is the proper frequency of the pulsating atom. 
The observer located at the receiver, in turn, assigns a frequency v R = to the 
pulsating atom. Consequently, 



vr = v /l-AMGK (m 2 r E ) 
ve y/l - AMGK (m 2 r R ) 

w 1 + 2MG [K (m 2 r R ) - K (m 2 r E )] . 

From this we immediately get the fractional shift 

— = VR ~ VE « 2MG [K (m 2 r R ) - K (m 2 r E )} . 
v ve 

Note that since K Q (x) is a monotonically decreasing function in the range < x < 
oo, is positive if r E > r R , and negative if r E < r R . Consequently, if the emitter is 
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nearer to the massive object than the receiver is, the shift is towards the red, but if the 
receiver is nearer the massive object, it is towards the blue. 

From the preceding considerations we come to the conclusion that the gravitational 
spectral shift is indeed a classical test of the BHT model. It can also be viewed, like in 
4D general relativity, as a direct test of the curvature of the spacetime. 




Emitter at 
fixed (re, Be) 



Receiver at 
fixed (r R ,6 R ) 



Figure 1. Spacetime diagram illustrating the worldlines of two successive 
identical signals. 



4. Gravitational time delay 

Another interesting effect that can be obtained from the linear approximation of new 
massive gravity is the time delay suffered by a light signal sent by an observer — situated 
at a fixed point in space in the gravitational field generated by a massive object — to 
a small object and reflected back to the observer. The small object is supposed to be 
located directly between the observer and the huge body (see Fig. 2). Consider, in this 
spirit, a light pulse that moves along a straight line connecting the observer and the 
small object. It is easy to show that the coordinate time for the whole trip (observer — > 
small object — > observer) is given by 



At 



G 




1 + 4MG[K (m 2 r) + 2 In 



1 - 4MGK (m 2 r) 



dr. 



(2* 



Accordingly, the proper time lapse measured by the observer, whose clock, of course, 
records proper time, has the form 
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Ar G = 2^1 - 4MGK (m 2 r 2 ) £ ^ 



1 + AMG 



K (m 2 r) + 21n^ 



AMGK (m 2 r) 



■dr. (29) 



0) 



Massive object 
of mass M 



Small object 
of mass m 



Observer 



Figure 2. Time delay in "radar sounding". 
On the other hand, the distance traveled by the light pulse is equal to 



2 f 2 Jl + AMG[K (m 2 r) + 2 In —}alr. 

Jri V r 

Consequently, on the basis of the classical theory we should expect a round-trip time of 



At, 



c 



1 + AMG 



K (m 2 r) + 2 In • 



dr. 



(30) 



From (29) and (30), we arrive to the conclusion that Atq ^ Arc. Note that in the 



m 2 — > oo limit, Arc — Arc — 2 y 1 + 8MGln ^dr, which clearly shows that there 
is no time delay in the framework of 3D general relativity, as expected. 
On the other hand, Eqs. (29) and (30), tell us that 



At, 
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V2 



1 + AMG (^K (m 2 r) + In 
- 4MG[K (m 2 r 2 )}(r 2 - ri ), 



dr 
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At g - At c « AMG 



K (m 2 r)dr - (r 2 - r 1 )K (m 2 r 2 ) 

n 

AMG [K (m 2 r ) - K (m 2 r 2 )] (r 2 - n), 



where ri < ro < r 2 . Hence, we come to the conclusion that there is a new- massive- 
gravity- induced time delay 

Ar G - Ar c w 4MG [if (m 2 r ) - ^ ("i 2 r 2 )] (r 2 - n). (31) 
5. Final remarks 

As is well-known, three-dimensional Einstein gravity without sources is physically 
vacuous because Einstein and Riemann tensors are equivalent in D = 3. In addition, 
the quantization of the gravity field does not give rise to propagating gravitons since the 
spacetime metric is locally determined by the sources. Consequently, the description 
of gravitational phenomena via 3D gravity leads to some bizarre results, such as the 
following. 

• Lack of a gravity force in the nonrelativistic limit. 

• Gravitational deflection independent of the impact parameter. 

• Complete absence of gravitational time dilation. 

• No time delay. 

It can be shown that the first two odd phenomena in the above list do not take place 
in the context of the BHT model [23]. In fact, in the framework of the latter, short-range 
gravitational forces are exerted on slowly moving particles; besides, the light bending 
depends on the impact parameter, as it should. On the other hand, the remaining 
strange phenomena in the aforementioned list, as we have shown, do not occur in the 
BHT system either. Indeed, both time delay and spectral shift do take place in the 
context of the new massive gravity. Like in 4D general relativity, gravitational time 
dilation and gravitational time delay are also tests of the BHT model. It is worth 
noticing that the basis for these tests is the time-independent solution of the linearized 
BHT field equations produced by a static spherical mass. 

One of the main reasons for studying 3D gravity models is in reality to try to find 
out a gravity system with less austere ultraviolet divergences in perturbation theory. 
Since general relativity in 3D is dynamically trivial, the BHT model, which is tree- 
level unitary, is an important step in this direction. This kind of research conducted 
in lower dimensions certainly helps us to gain insight into difficult conceptional issues, 
which are present and more opaque in the physical (3+l)-dimensional world. Another 
strong argument in favor of considering massive gravity theories, as we have already 
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commented, is the fact that the present accelerated expansion of the universe could be 
partially attributed to a graviton mass-like effect. 

It is worth mentioning that the triviality of 3D general relativity can also be cured 
by adding to the EH action in 3D a parity-violating Chern-Simons term. The resulting 
model is usually known as topological massive gravity (TMG) [211 [25] . Nonetheless, 
in contrast with TMG, 3D massive gravity has the great advantage of being a parity- 
preserving theory. On the other hand, since 3D higher-derivative gravity (3DHDG) - 
which is defined by the Lagrangian density £3dhdg = \f§ (ff R + § Rp U + f -^ 2 ) — i s 
nonunitary at the tree level [26], it would be interesting to verify whether the addition 
of a topological Chern-Simons term (£ C s = § eAMI T p CTA [d M r CT + fr" 7 r w ], where fi is 
an arbitrary parameter) to this higher-order model would cure the nonunitarity of the 
former. It can be shown that in order to avoid ghosts and tachyons in the mixed theory 
(£ = £3dhdg + £cs) ^ ne f ouowrn g constraints on the parameters must holdjj] [27] : 

(spin-2 sector) : a < 0, (3 > 0, 
(spin-0 sector) : a > 0, 3/3 + 8a > 0. 

Therefore, for arbitrary values of the parameters, the model at hand is nonunitary 
at the tree level, which clearly shows that the topological Chen-Simons term is not a 
panacea for 3DHDG's unitarity problem. Nevertheless, if we prevent the spin-0 mode 
from propagating by choosing 3(3 + 8a = 0, the resulting model is tree-level unitary. 
It is amazing that the above condition is exactly the same constraint that appears 
in the BHT model (m — > oo limit). We call attention to the fact that, contrary to 
popular belief, the addition of a Chern-Simons term to a tree-level unitary model is 
not necessarily a guarantee that the resulting model will be tree-level unitary [26]. For 
instance, the addition of a Chern-Simons term (Ccs) to three-dimensional R + aR 2 
gravity (C^+aR 2 = ( — + ^f~)\/#)> which is tree-level unitary, spoils the unitary of 
the latter [26]. Therefore, in some cases the coexistence between the topological Chern- 
Simons term and 3D higher-derivative gravity theories is conflicting. 

To conclude we remark that recently the nonlinear classical dynamics of the BHT 
model was exhaustively investigated by de Rham, Gabadadze, Pirtskhalava, Tolley 
and Yavin [28], who found that the theory passed remarkably nontrivial checks at the 
nonlinear level, such as the following. 

• In the decoupling limit of the theory, the interactions of the helicity-0 modes are 
described by a single cubic term, the so-called cubic Galileon [29]. 

• The conformal mode of the metric coincides with the helicity-0 mode in the 
decoupling limit. 

• The full theory does not lead to any extra degrees of freedom, which suggests that 
a 3D analog of the 4D Boulware-Deser ghost is not present in the BHT system. 

| The massless excitation, like the massless excitation of 3D general relativity, is a not a dynamical 
degree of freedom, i.e., it is nonpropagating. 
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Appendix A. Propagator 

In order to find the propagator related to the Lagrangian density in Eq. (1) it is very 
convenient to work in terms of the Barnes-Rivers operators in the space of symmetric 
rank-two tensors. The complete set of 3-dimensional operators in momentum space is 
[301131] 



(A - 2) 

P au,K\ = ^H^kX, (A. 3) 



fiu,nX 
,(0-sm) 



1 



p*::r = ( A - 4 ) 



where 8^ u = 77^ — and = -p^ are, respectively, the usual transverse and 
longitudinal projection operators. The multiplicative table for these operators is 
displayed in Table I. 
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Table Al. Multiplicative table for the Barnes-Rivers operators 
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p(0— su>) 





p(0-s) 


p(0—ws) 








p(0—ws) 





p(Q-to) 






To compute the graviton propagator we need the bilinear part of the Lagrangian 
density (1). With the gauge fixing -^(d^ u ) 2 (de Donder gauge), and going over to 
momentum space we reproduce (5). The task of computing the operator O is greatly 
facilitated if we appeal to the following identities 
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[p(2) + pd) + p(0-.) + p(0-»)] ^ = 1 ( ^ A + 



2p(0-s) _|_ p(0-w) _|_ v /2(p(°~ slu ) _|_ p(O-ws)^ 



[2P (1) + 4P (0 -"' ) ] = ^(rj^k u k x + T^fc* + ^a&A + V^k x ), 



2p(o-w) _|_ v /2(P ( - ~ sw - ) + p(°- ws ) 



Now, if we write the operator in the generic form 

= Xl P {1) + x 2 P (2) + x s P { °- s) + x w P {0 ' w) + x S5i ,P (0 - su,) + x ws P ( °- ws \ 
and take into account that OO^ 1 = I, where O^ 1 is the propagator, we promptly find 



0-1 = lp(D + lp(2) + _L 

Xl X2 XgXyj X SW X WS 



T p(O-s) 1 p(O-w) 



_ r p(O-sw) _ p(0— u>s)l 

From (A. 6) and (5) we obtain (6). 



(A.6) 



Appendix B. A useful result 

Theorem 1 If m is the mass of a generic physical particle related to a given 3D 
gravitational model and k is the corresponding exchanged momentum, then 

(T^-V)| fe2=m2 >0 and (T 2 „-T 2 )| fc2=0 = 0. 

Here T^ v [= T u ^) is the external conserved current. 

We begin by remarking that the set of independent vectors in momentum space, 
k^ = (&°,k), k^ = (A; ,— k), e = (0, e), where e is a unit vector orthogonal to k, is a 
suitable basis for expanding any three- vector V^(k). Using this basis we can write the 
symmetric current tensor as follows 

T i* = AWk v + B ~ k »~ k v + Ce » e v + Dk (p~ k u) + Ek^e v) + Fk^e v \ 

where aH u) = \(a»b v + Va"). 

The current conservations gives the following constraints on the coefficients A, B, 
C, D, E, and F: 



Ak 2 + ^(k 2 + k 2 )=0 

B{kl + k 2 ) + jk 2 = 
Ek 2 + F(k 2 + k 2 ) =0 



(B.l) 

(B.2) 
(B.3) 
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From Eqs. (Bl) and (B2), we get Ak A + B(k 2 + k 2 ) 2 , while Eq. (B3) implies 
E 2 > F 2 . On the other hand, saturating the indices of T^ v with momenta we arrive 
at a consistent relation for the coefficients A, B, and D: 

Ak 4 + B(k 2 + k 2 ) 2 + Dk 2 (k 2 + k 2 ) = 0. (B.4) 

After a lengthy but otherwise straightforward calculation using the earlier 
equations, we obtain 



rp2 rp2 



rp2 rpl 

1 jJLV 



Therefore, 



{T 2 



k 2 (A-B) _ cr 2 

V2 V2 



k 



+ 'j(E 2 -F 2 ), 



-(E 2 - F 2 ) — 2C(A - B) 



(B.5) 



l -T 2 )\ k 2 =m > > Q and {T 2 V — T 2 )| fc 2 =0 = 0. 
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